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Abstract. 

We present a new set of multi-million particle SPH simulations of the for- 
mation of disk dominated galaxies in a cosmological context. Some of these 
galaxies are higher resolution versions of the models already described in Gover- 
nato et al (2007). To correctly compare simulations with observations we create 
artificial images of our simulations and from them measure photometric Bulge 
to Disk (B/D) ratios and disk scale lengths. We show how feedback and high 
force and mass resolution are necessary ingredients to form galaxies that have 
flatter rotation curves, larger I band disk scale lengths and smaller B/D ratios. 
A new simulated disk galaxy has an I-band disk scale length of 9.2 kpc and a 
B/D flux ratio of 0.64 (face on, dust reddened). 



1. The Simulation Campaign 

Here we present results from an ongoing project aimed at simulating a repre- 
sentative sample of gal axies that forin in a ACDM cosmology. Simulations are 
run with G ASOLINE (IW adsley etaLl l2004l ) and use the "blastwave feedback" 



described in IStinson et al.i (,2006 ). Three of the four models presented here use 
the same initial conditions as in Governato et al. (2007, G07) and thus the 
galaxy halos have the same merging history and form inside the same large scale 
structure. However, these new simulations use 8 times more resolution elements 
and a 50% smaller force softening than those presented in G07. To understand 
the numerical limitations of our work and to better compare with the existing 
literature, we vary the number of resolution elements and the feedback imple- 
mentations. In one case (named MWl) we run three models with the number of 
resolution elements varying from a few tens of thousand to a few million particles 
and force resolution from 1.2kpc to 0.3kpc. To allow a straightforward compari- 
son between simulations and observations, the galaxies' B/D ratios and disk scale 
lengths are measured using GALFIT on d ust reddened, face-on I and K band 
artificial images produced using SUNRISE (|Jonssonll2006l ). Recent works by dif- 



ferent groups and using different hydrodynamical codes have highlighted the role 
of feedback in creating gala.xies with a dominant stellar component supported 
by rotation ( Okamoto et al.ll2005l : [Robertson et al.l 20041 ). Significant improve- 



ments have been made in reproducing a number of scaling properties of late type 
galaxies: namely the age - circular velocity, the Tully Fisher, the stellar mass 



metallicity relations and the abundance of ga laxy satellites (| Governato et al 



20071 : iBrooks et aDl2007l : iLibeskind et~alll2007l ). However, a well known prob 



lem of most galaxies formed in cosmological simulations is a massive, centrally 
concentrated spheroidal component that creates unrealistically declining rota- 
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tion curves ( Eke et alll200ll : [Governato et al.ll2007l ). The explanations proposed 
in the hterature invoke both physical and numerical arguments: catastrophic 
angular momentum loss due to dynamical friction suffered by gas rich subhalos 



(jNavarro et al.lll996l '). spurious gravitational torques between the disk and the 



surrounding halo of dark matter and gas caused by the noisy mass distribution 



at low resolution ( Kaufmann et al. 20071). or artificial pressure gradients at the 



cold/hot gas interface ( Okamoto et alj 20031 ). 




Figure 1. Left: The kinematic disk, bulge and halo stellar mass fraction 
(blue, red and orange) in the MWl-HR, MWl-MR and MWl-MR-TH models. 
All runs use the "blastwave" feedback, but the one marked "TH", where 
"thermal" feedback is instead used. MR is the close equivalent of the MWf 
galaxy described in G07. Histograms are normalized to 100%. The total 
stellar masses differ by less than fO%. Right: The circular velocity profile 
from the total mass distribution (Vc = (M/r)^/^) for the same galaxy halo, 
run with different resolutions and feedback implementations. The circular dot 
marks the radius that contains 85% of the stars. HR, MR and LR runs use a 
total of 5 million, 600k and 70k DM, gas and star particles within the virial 
radius. 



Run I-band (kpc) Kinem. B/D u -r 



MWf-HR 


3.5 


0.83 (0.38) 


1.6 


MWl-MR 


3.3 


1.56 


1.7 


MWl-MR-TH 


2.1 


1.51 


1.6 


MWl-LR 


2.1 


1.35 


2.0 



Table 1. Summary of the photometric properties of the MWl run obtained 
using SUNRISE and GALFIT. The value in brackets is the photometric B- 
band B/D ratio for MWl-HR. MWl has a total mass of lO^^ M© 



SN feedback is a promising solution as it decoup l es ba ryons from the dark 

( 20071 ) recently presented 



Naab et al. 



matter (Zavala et al., 2007). However, 
results from simulations that did not include feedback from SNe: the resulting 
circular velocity profile of their simulated galaxies changed significantly with 
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Figure 2. Left: The rotation curves for the higher resolution version equiv- 
alent of the models in G07. The dashed line is the curve for MWl published 
in the paper, equivalent to the new "MR" model. The dotted line is the low 
resolution MWl-LR run. The vertical lines mark the radius at which Yrot 
is measured: 3.5 Tband R^s. This is a radius were the amount of included 
mass (DM, stars & gas) is close to converging (Fig. lb) and gives more reli- 
able results compared to measuring V^ot at 2.2 R^ where numerical effects 
might still be present. Right: the TuUy Fisher relation for a sample of four 
high resolution runs vs. data in Giovanelli et al 97) Large filled dots: SSPs 
from Girardi et al (2000). Smaller filled dots: SSPs from Starburst99. Open 
square: MWl-LR. 



increasing resolution, with the peak velocity dropping from 320 km/sec to 220 
when the numb e r of p articles was increased from 40^ to 200'^ (their Fig.l). 
iKaufmann et al.l ( 20071 ) used simulations of isolated galaxies to show that cold, 
rotationally supported disks lose a large fraction of thier initial angular mo- 
mentum if the number of gas particles is below 10®. The literature presents 
encouraging evidence that both a physically motivated description of feedback 
and a high number of resolution elements help to form galaxies with at least 
some of the properties of the real ones. However, an analysis of a uniform sam- 
ple of simulations with techniques directly comparable to those used by observers 
(who measure the light rather than the underlying mass distribution) has been 
lacking, making a comparison between observations and theoretical models more 
open to interpretation. 



2. Results from Simulations 

To show the separate effects of feedback and resolution we focus on a simulation 
of a Milky way like galaxy (total halo mass 10^^ M0). The disk, halo and 
bulge components are first identified based on their kinematics (Brook et al 08, 
in prep). All runs adopted the "blastwave" feedback, but those n amed "TH", 
that use the much less effective "thermal" feedback ( Katd Il992l ) which is a 
lower limit to the possible effects of feedback (the two methods inject the same 
amount of energy into the ICM). All runs include the effects of a uniform cosmic 
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Figure 3. Sunrise Unreddened I-band images of the MWl model. Up- 
per panel: High to low resolution (with blastwave feedback). Lower panel: 
The medium resolution model with blastwave feedback, no SN feedback and 
"Thermal" feedback. Galaxies with feedback and high resolution form disks 
with larger scale lengths and/or less concentrated bulges. The dust reddened, 
photometrically decomposed B/D ratio in B and I bands is lower than the 
one measured using the kinematics. 



UV background. The kinematically defined stellar disk component becomes 
more dominant as resolution increases or a physically realistic description of 
feedback effects is adopted (Fig. la). Bulge masses decrease significantly in the 
high resolution run with feedback (HR vs MR and & MR-TH). Feedback is also 
particularly effective in reducing the fraction of halo stars, a consequence of a 
smaller stellar mass in stripped galaxy satellites (MR vs. MR-TH) (Brook et al 
04). Fig. lb shows the circular velocity profile (Vc(r) = yjM/R) for the same 
galaxy as resolution and feedback are varied and shows results consistent with 
those in lNaab et al. ( 20071 ) and Kaufmann et al (2007): higher resolution models 
are less centrally concentrated. While this plot emphasizes the role of resolution, 
we have verified that effect of feedback on the mass distribution is larger in less 
massive halos with a shallower potential well, consistent with Zavala et al (2007). 
After a bulge/disk decomposition with GALFIT, photometric, dust reddened I- 
band disks have larger scale lengths (Rd) and galaxies have bluer colors in higher 
resolution runs that include the blastwave feedback (Table 1 and Fig. 3). Fig. 2a 
shows the rotation curves from cold gas or young stars for the high resolution 
version of the galaxies in G07. Compared with the lower resolution models in 
G07 rotation curves are flatter, as the very central regions are less dense and the 
disks dynamically colder in the outer parts. However, they are still not as flat 
as those of most real galaxies of similar mass. Fig. 2b shows the resulting I-band 
Tully Fisher relation. Yrot is measured at 3.5 R^- There is good agreement 
between data and the high resolution models, confirming results in G07, while 
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the low res MWl model is significantly shifted to the right compared to its high 
resolution equivalent and the observed distribution. We added to the sample a 
recently simulated large spiral galaxy with an R^^ of 9.2, 9.2 and 9.1 kpc in the 
B,I and K bands respectively. Its dust reddened I-band B/D is = 0.64. 
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